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Abstract 

The ball (or sphere) packing problem with equal balls, without any sym- 
metry assumption, in a 3-dimensional space of constant curvature was set- 
tled by Boroczky and Florian for the hyperbolic space in [2] and by 
proving the famous Kepler conjecture by Hales 1 8] for the Euclidean space 
E^. The goal of this paper is to extend the problem of finding the densest 
geodesic ball (or sphere) packing for the other 3-dimensional homogeneous 
geometries (Thurston geometries) 

S^xR, H^xR, SL^, Nil, Sol, 

where a transitive symmetry group of the ball packing is assumed, one of 
the discrete isometry groups of the considered space. 

*AMS Classification 2010: 52C17, 52C22, 53A35, 51M20 
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Moreover, we describe a candidate of the densest geodesic ball packing. 
The greatest density until now is ^ 0.85327613 that is not realized by pack- 
ing with equal balls of the hyperbolic space H^. However, it attains e.g. at 

3 

horoball packing of H where the ideal centres of horoballs lie on the ab- 
solute figure of H'^ inducing the regular ideal simplex tiling (3, 3, 6) by its 
Coxeter-Schlafli symbol. In this work we present a geodesic ball packing in 
the X R geometry whose density is 0.87499429. The extremal config- 
uration is described in Theorem 2.8, Our conjecture and further remarks are 
summarized in Section 3. 



1 Preliminary results 

1.1 Geodesic ball packings in spaces of constant curvature 

Finding the densest (not necessarily periodic) packing of balls in the 3 -dimen- 
sional Euclidean space is known as the Kepler Problem: No packing of spheres 
of the same radius has a density greater than the face-centered cubic packing. 
This density can be realized by hexagonal layers (in continuum many ways). 
This conjecture was first published by Johannes Kepler in his monograph The 
Six-Cornered Snowflake (1611), this treatise inspired by his correspondence with 
Thomas Harriot (see Cannonball Problem). In 1953, Laszlo Fejes Toth reduced 
the Kepler conjecture to an enormous calculation procedure that involved spe- 
cific cases, and later suggested that computers might be helpful for solving the 
problem. In this way the above four hundred year mathematical problem has fi- 
nally been solved by Thomas Hales |[8l. He had proved that the guess of Kepler 
from 1611 was correct. In mathematics sphere packing problems concern arrange- 
ments of non-overlapping equal spheres (rather balls) which fill a space. Space 
is the usual three-dimensional Euclidean space. However, ball (sphere) packing 
problems can be generalized to the other 3-dimensional Thurston geometries, but 
a difficult problem is - similarly to the hyperbolic space - the exact definition of 
the packing density. 

In an n-dimensional space of constant curvature E", H", S*^ (n > 2) let dn{r) 
be the density of n + 1 spheres of radius r mutually touching one another with 
respect to the simplex spanned by the centres of the spheres. L. Fejes Toth and 
H. S. M. Coxeter conjectured that in an n-dimensional space of constant curvature 
the density of packing balls of radius r can not exceed dn{r). This conjecture has 
been proved by C. Rogers in the Euclidean space E*^ fTSl. The 2-dimensional 
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spherical case was settled by L. Fejes Toth in [|71 and in [[T| K. Boroczky proved 
the following generalization: 

Theorem 1.1 (K. Boroczky) In an n-dimensional space of constant curvature 
consider a packing of spheres of radius r. In spherical space suppose, that r < |. 
Then the density of each sphere in its Dirichlet-Voronoi cell cannot exceed the 
density ofn + 1 spheres of radius r mutually touching one another with respect to 
the simplex spanned by their centres. 

Remark 1.2 

1. In the hyperbolic space this result can be extended to r = oo /[7]/ where the 
densest horoball packing can be realized by different regular arrangements ^Tl\l . 

2. If we allow horoballs of different types at the various vertices of a totally asymp- 
totic simplex and generalize the notion of the simplicial density function in H", 
(n > 2) then the Bdrdczky-Florian type density upper bound does not remain 
valid for the fully asymptotic simplices / |25l/ . / [2^ . 

1.2 Geodesic ball packings under discrete isometry groups 

E. Molnar has shown in [13], that the homogeneous 3-spaces have a unified in- 
terpretation in the projective 3-sphere VS^iY'^, V4, M). In our work we shall use 
this projective model of each Thurston geometry and the Cartesian homogeneous 
coordinate simplex Eo{eo),E^{ei),E^{e2), E^ie^), ({6^} C with the unit 
point E{e = eg + ei + 62 + 63)) which is distinguished by an origin Eq and by the 
ideal points of coordinate axes, respectively. Moreover, y = cx with < c G M 
(or c G R \ {0}) defines a point (x) = (y) of the projective 3-sphere VS^ (or that 
of the projective space where opposite rays (x) and (— x) are identified). The 
dual system {(e*)} C V4 describes the simplex planes, especially the plane at 
infinity (e°) = EfE'^Ef, and generally, v = defines a plane [u) = (v) of 
VS^ (or that of V^). Thus = x-u = yv defines the incidence of point (x) = (y) 
and plane (u) = (v), as (x)I('u) also denotes it. Thus the Thurston geometries 
can be visualized in the affine 3-space A'^ (so in E'^) as well. 
Let X be one of the remaining 5 Thurston geometries 

S^xR, H^xR, SlS., Nil, Sol 



where the geodesic curves are generally defined as having locally minimal arc 
length between their any two (near enough) points. The equation systems of the 
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parametrized geodesic curves j{t) in our model can be determined by the general 
theory of Riemann geometry. Then geodesic sphere and ball can be usually de- 
fined as follows below. We consider only those geodesic ball packings which are 
transitively generated by discrete groups of isometries of X and the density of the 
packing is related to its Dirichlet-Voronoi cells. 

Definition 1.3 The distance d{Pi, P2) between the points Pi ^ X and P2 & X is 
defined by the arc length of the geodesic curve from Pi to P2. 

Definition 1.4 The geodesic sphere of radius p (denoted by Spj^{p)) with centre 
at the point Pi is defined as the set of all points P2 in the space with the condi- 
tion d{Pi, P2) = p. Moreover, we require that the geodesic sphere is a simply 
connected surface without selfintersection of space X. 

Definition 1.5 The body of the geodesic sphere of centre Pi and of radius p in 
space X is called geodesic ball, denoted by Bp^{p), i.e. Q e Bp^{p) iff Q < 
d{Pi,Q)<p. 

In the following let F be a fixed group of isometries of X. We will denote by 
d{Pi, P2) the distance of two points Pi, P2 (see Definition (1.3)). 

Definition 1.6 We say that the point set 

V{K) = {P e X : d{K,P) < d(K^,P)forallge T} 
is the Dirichlet-Voronoi cell (D-V cell) to T around the kernel point K E X. 
Definition 1.7 We say that 

Tp = {geT : P^ = P} 
is the stabilizer subgroup ofPeX in F. 

Definition 1.8 Assume that the stabilizer Tk = I the identity, i.e. F acts simply 
transitively on the T-orbit of K G X. Then let Bk denote the greatest ball of 
centre K inside the D-V cell T>{K), moreover let p{K) denote the radius of Bk. 
It is easy to see that 

piK) = min -d{K,K^). 
ger\i 2 
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Definition 1.9 If the stabilizer Tk > I then F acts multiply transitively on the 
V-orbit ofK&X. Then the greatest ball radius o/Bk is 

piK) = mill -d(K, K^) 
ger\rK 2 

where K belongs toaO- 1- or 2 -dimensional region ofX (vertices, axes, reflection 
planes). 

In both cases the F-images of form a ball packing with centre points K^. 



Definition 1.10 The density of ball packing B^ is 

VoI{Bk) 



5{K) 



VolV{K) ' 



It is clear that the orbit and the ball packing have the same symmetry 
group, moreover this group contains the starting crystallographic group F: 

SymK^ = SymB]^ > ^■ 



Definition 1.11 We say that the orbit and the ball packing is characteris- 
tic if SymK^ — F, else the orbit is not characteristic. 

1.2.1 Simply transitive ball packings 

Let F be a fixed group of isometrics in space X. Our problem is to find a point 
K e X and the orbit for F such that F^ = I and the density 5{K) of 
the corresponding ball packing B^{K) is maximal. In this case the ball packing 
B^{K) is said to be optimal. 

Our aim is to determine the maximal radius p{K) of the balls, and the maximal 
density 6{K). The considered space groups could have free parameters, so we 
have to find the densest ball packing for fixed parameters piV), then we have to 
vary them to get the optimal ball packing 

5{V) = max (S(K)). (1.1) 
K, p(r) 

We look for the optimal kernel point in a 3 -dimensional region, inside of a funda- 
mental domain of F. 
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1.2.2 Multiply transitive ball packings 

Similarly to the simply transitive case we have to find a kernel point K G X and 
the orbit for T such that the density S{K) of the corresponding ball packing 
B^{K) is maximal but here 7^ I- This ball packing B^{K) is called optimal, 
too. In this multiply transitive case we look for the optimal kernel point K in pos- 
sible 0- 1- or 2-dimensional regions C, respectively. Our aim is to deteremine the 
maximal radius p{K) of the balls, and the maximal density S{K). The considered 
space group can have also free parameters p{T), then we have to find the densest 
ball packing for fixed parameters, and vary them to get the optimal ball packing. 

6(T) = max (6(K)) (1.2) 
K&c,p{r) 

1.2.3 Packings in Nil space 

W. Heisenberg's famous real matrix group provides a non-commutative transla- 
tion group of an affine 3-space. The Nil geometry, which is one of the eight 
homogeneous Thurston 3-geometries, can be derived from this matrix group fT3\. 

In [I2TII I have investigated the geodesic balls of the Nil space and computed 
their volume, introduced the notion of the Nil lattice. Nil parallelepiped and the 
density of the lattice-like ball packing. Moreover, I have determined the densest 
lattice-like geodesic ball packing by a type of Nil lattices. The density of this 
packing is ^ 0.78085, may be surprising enough in comparison with the analo- 
gous Euclidean result ^ 0.74048. The kissing number of every ball in this 
packing is 14. By my conjecture the densest geodesic ball packing belongs to the 
above ball arrangement in Nil space. The symmetry group of this packing has 
also been described in [[T4|. 

1.2.4 Packings in x R space 

This Seifert fibre space is derived by the direct product of the hyperbolic plane 
and the real line R. In ETll I have determined the geodesic balls of x R space 
and computed their volume, defined the notion of the geodesic ball packing and its 
density. Moreover, I have developed a procedure to determine the density of the 
simply or multiply transitive geodesic ball packings for generalized Coxeter space 
groups of X R and apply this algorithm to them. For the above space groups 
the Dirichlet-Voronoi cells are "prisms" in x R sense. The optimal packing 
density of the generalized Coxeter space groups is: ^ 0.60726. I'am sure, that in 
this space there are denser ball packings. 
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Remark 1.12 1. So far there are no results for the geodesic ball packings in 
Sol and SL2R geometries (by my knowledge). 

2. In Nil and Sol spaces I have studied the so-called translation ball packings 
(see HITj . ^22\1) but I do not consider these cases in this work. 

2 On X R space 

The S^ X R geometry can be derived by the direct product of the spherical plane 
S^ and the real line R. In [4] J. Z. Farkas has classified and given the complete 
list of the space groups in S^ x R. The S^ x R manifolds up to similarity and 
diffeomorphism were classified by E. Molnar and J. Z. Farkas in llH. In [|23|| 
I have studied the geodesic balls and their volumes in S^ x R space, moreover I 
have introduced the notion of geodesic ball packing and its density and determined 
the densest simply and multiply transitive geodesic ball packings for generalized 
Coxeter space groups of S^ x R, respectively. The density of the densest packing 
is ^ 0.82445. 

In paper ll24l I have studied the simply transitive locally optimal ball pack- 
ings to the S^ X R space groups having Coxeter point groups and at least one of 
the generators is a non-trivial glide reflection. I have determined the densest sim- 
ply transitive geodesic ball arrangements for the above space groups, moreover I 
computed their optimal densities and radii. The density of the densest packing is 
^ 0.80408. 

Now, we shall discuss the simply and multiply transitive ball packings to a 
given space group. But let us start first with the necessary concepts. The points 
in S^ X R geometry are described by (P, p) where P E and p G R. The 
isometry group /som(S^ x R) of S^ x R can be derived by the direct product of 
the isometry group of the spherical plane /som(S^) and the isometry group of the 
real line /som(R). 

The structure of an isometry group T C /som(S^ x R) is the following: T := 
{{Ai X pi),...(yl„ X where Ai x pi := Ai x (Pi,^) := {gi,ri), (i G 

{1, 2, . . . n} and A^ G /som(S^), Ri is either the identity map 1r of R or the 
point reflection 1r. gi := A^ x is called the linear part of the transformation 
{Ai X Pi) and is its translation part. The multiplication formula is the following: 



(A X R^,n) o (A2 X R2,r2) = {{A^A2 x R.R^^nR^ + r^). (2.1) 
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Definition 2.1 A group of isometries V C /som(S^ x R) is called space group 
if the linear parts form a finite group Vq called the point group of V, moreover, 
the translation parts to the identity of this point group are required to form a one 
dimensional lattice Lr o/R. 

Remark 2.2 It can be proved, that the space group V exactly described above has 
a compact fundamental domain J-p- 

We characterize the spherical plane groups by the Macbeath-signature (see fT2], 

m). 

In this paper we deal with a class of the x R space groups 4q. I. 2 (with a 
natural parameter q > 2, see /|?|/ ). Each of them belongs to the glide reflection 
groups, i.e. the generators g,, (i = 1, 2, . . . m) of its point group Tq are reflections 
and at least one of the possible translation parts of the above generators differs 
from zero ( see Il24\l ). 

2.1 Geodesic curves and balls in x R space 

In ll23l and [24] I have described the equation system of the geodesic curve and 
so the geodesic sphere: 

x{t) = e^''™''cos (rcos-u), 
?/(r) = e^^^^^ sin (r cos v) cos u, 
z[T) = e sm (r cost;) smw, ^ 

TT TT 
— TT < M < TT, < ti < — 

- ' 2 - - 2 

of radius p = r > 0, centre x(0) = 1, y{0) = 0, z{0) = 0, and with longitude u, 
altitude v, as geographical coordinates. 

In f23\ I have proved that geodesic sphere S{p) in x R space is a simply 
connected surface in if and only if p G [0, vr), because if p > tt then there is 
at least one v G [— |, |] so that ?/(r, v) = z{t, v) = 0, i.e. selfintersection would 
occur. Thus we obtain the following 

Proposition 2.3 The geodesic sphere and ball of radius p exists in the x R 
space if and only if p E [0, tt). 

We have obtained (see [|23l ') the volume formula of the geodesic ball B{p) of ra- 
dius p by the metric tensor gij and by the Jacobian of (2.2) and a careful numerical 
Maple computation for given p by the following integral: 
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Theorem 2.4 



voi{B{p)) = I — ^ 2\3/2 ^^ 



2 



|r ■ sinfcosfvlr)! dw dr 



2tv I I \t ■ sm(cos(v)T)\ dv dr. 







2 



(2.3) 



The fundamental domain of the studied space groups can be combined as a 
fundamental domain of the spherical group with a part of a real line segment. 
This domain is called x R prism. In [|23l we have shown the following 



Q 




Co 



Figure 1 : Trigonal and digonal prisms 



Theorem 2.5 The volume of a trigonal prism VbqBiB2CqCiC2 <^f^d of a digonal 
prism VboBiCqC'i in x R (see Fig. l.a-b) can be computed by the following 
formula: 

VoliV) =A-h (2.4) 



where A is the area of the spherical triangle AqAiA2 or spherical digon AqAi in 
the base plane 11 belonging to the fibre coordinate t = and h = BqCq is the 
height of the prism or digon, respectively. 
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2.2 Optimal ball packings for the space group 4q. 1. 2, (q=2) 

We consider a x R space group (see flU |23l |24l| with point group Fq generated 
by three reflections gj, = 1, 2, 3) 

(+, 0, [ ]{i2,2,q)}), q>2, 

To = (gl, g2, g3 - g?, gl gl, (glg3)^ (g2g3)^), (glg2)^)- 

The possible translation parts ti, T2, of the corresponding generators of Tq 
will be determined by (2.1) and by the defining relations of the above point group. 
Finally, we obtain six non-equivariant solutions from the so-called Frobenius con- 
gruence relations: 

(ri,r2,r3) = (0,0,0), (0,0,^), (^44^' ^^4'°^' ^^'I'l^' 

If (ti, T2, ts) = (0, 0, |) then we have obtained the x R space group 4q. L 2 
(for a fixed q,2 < q eN). 

The fundamental domain of the point group of the considered space group is 
a spherical triangle with angle ^, |, | in the base plane H. It can be 

assumed, that the fibre coordinate of the center of the optimal ball is zero and it is 
an interior point of A1A2A3 triangle (see Fig. 2). 

In the following we consider ball packings, only to q = 2. We shall apply 
the above introduced Cartesian homogeneous coordinate system (see Fig. 2) and 
the usual geographical coordinates {(j),9), (— vr < (p < n, — | < 6* < |) of the 
sphere with the fibre coordinate t G R. 

2.2.1 The simply transitive ball packing 

We consider an arbitrary interior point K{x'^, x^,x'^, x^) = K{(j), 6) of spherical 
triangle AiA2A^ in the above coordinate system in our model (see Fig. 2). Our 
aim is to determine the maximal radius R of the balls, and the maximal density 
S{K, R) where K is an inner point of A1A2A3 triangle or K E A1A2. 

Thus, q = 2, = T, the x R group is 42.1.2 now, the R-translation is 2r. 
(Instead of the above 6 cases we have only 4 space groups for the corresponding 
point group.) The fundamental domain of the point group of the considered space 
group is a spherical triangle A1A2A3 with angles |, |, | in the base plane 11 (see 
Fig. 2). It can be assumed by the homogeneity of x R, that the fibre coordinate 
of the center of the optimal ball is zero. 

x^ = 1, x^ = cos (f) cos 6 , x'^ = sin (j) cos 9 , x^ = sin^ (2.5) 
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Let Br{R) denote a geodesic ball packing of x R with balls Bk{R) of 
radius R (to be determined yet) where their centres K give rise to the orbit . 
In the following we consider to each ball packing the possible smallest translation 
part r(K, R) (see Fig. 2) depending on F, K and R. A fundamental domain of F 
Tr{K,R) is the corresponding Dirichlet-Voronoi cell where its volume is equal to 
the volume of the prism which is given by the fundamental domain of the point 
group Fo of F and by the height 2t{K, R). (It is clear that the optimal ball Bk has 
to touch some faces of its D-V cell to F around the kernel point K.) The images 
of J^t{k,r) by discrete isometry group F covers the x R space without overlap. 
For the density of the packing it is sufficient to relate the volume of the optimal 
ball to that of the solid ^t{k,r)- Analogously to the Euclidean case the density 
S{R, K) of the geodesic ball packing Br{R) can be defined : 

Definition 2.6 

^.MftWll^. (2.6) 

Vol{J'riK,R)) 
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Figure 3: Balls of the simply transitive optimal ball arrangement with base plane 
TT which is represented by the spherical triangle ^41^2^3 in figure b. 



(2-V) 



The optimal ball arrangement Bopt{K, R) has to satisfy the following equations: 

(a) d{K,K'^') = 2R= d{K,K3^), 
{b) d{K, K^^') = 2R< d{K, K^asrgs-^ < AR = d{K, K'^) 

We consider two main ball arrengements: 

1. We denote by i3r(^i, Ki) those packing where requirements (2.7) and 
d{K, K^^) = 2R hold (see Fig. 2). 

2. We denote by By{R2, K2) those packing where requirements (2.7) and 

d{K, K^aargs^ = (i{K, K^^) = AR hold (see Fig. 2). 

First we determine the coordinates of the points K^, (i = 1, 2) (K^ with parame- 
ters (pi and 6i), the radius Ri of the ball, the volume of a ball B{Ri) and the density 
of the packing in both main cases. We get the following solutions by systematic 
approximation, where the computations were carried out by Maple V Release 10 
up to 30 decimals: 

01 = - ^ 0.78539816, 9^ ^ 0.55737781. Ri ^ 0.64360446, 

4 ' ' (2.8) 

Vol{B{Ri) ^ 1.08624788, d{Ri,Ki) ^ 0.53722971. 
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Figure 4: Some neighbouring balls of the simply transitive optimal ball arrange- 
ment R)=BriRi,Ki) 

02 = - ^ 0.78539816, ^2 = 0, R2 = - ^ 0.78539816, 

4 ' ' 4 (2.9) 

Vol{B{R2) ^ 1.94735865, 5(/?2,i^2) ~ 0.39461737. 

We obtain by careful discussion of the density function 6{R, K) (R E [Ri, R2]) 
of the considered ball packing the following: 

Theorem 2.7 The ball arrangement Br{Ri, Ki) (see Fig. 3.a-b and Fig. 4) pro- 
vides the densest symply transitive ball packing belonging to the x R space 
group 4q. I. 2 (g = 2). 

2.2.2 The multiply transitive ball packings 

To determine the optimal multiply ball packing we have to study 2-cases: 

1. K is an inner point of the spherical geodesic segment (or AiA^i). 

In this situation the point K and its images by F = 4q. \.2 [q = 2), as 
the centers of the optimal ball arrangement Bopt{K,R) have to hold the 
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following requirements because of an arbitrary ball of the optimal packing 
is fixed by its neighbouring balls: 

(a) d{K,K^-)=2R = d{K,K^^-^), 
(b) 2R < d{K, K^a^rg^-j = ^(/^^ x'^r-^ < 4^, ^ " 

It is easy to see that in this case we get the optimal packing if K = A2 (or 
K = Ai) with the following data: 

(j)^ = -^ 1.57079633, 6*3 = 0, R3 = - ^ 1.57079633, 

2 ,3,^2 (2.11) 

Vol{B{Rs)) ^ 13.74539472, 5(i?3,i^3) ~ 0.69634983. 

K = A^. The Fig. 5 shows the orbit of the point ii' = by the considered 
space group. The images of K lie on a line through the origin and ^43. 
Fig. 6 and Fig. 7. show this interesting ball arrangement whose data are the 




Figure 5: The orbit of = A3 by the group T = 4q. I. 2 (g = 2). 
following 
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04 = _ ^ 0.78539816, ^4 = 2 ~ 1-57079633, ~ 1.81379936, 
Vol{B{Ri)) ^ 20.00238509, d{Ri,Ki) ^ 0.87757183. 

(2.12) 

The "outwards transformed" images of a balls sourround the previous balls 
(see Fig. 7.) thus the touching number of this packing is 4. Finally, we get 
the following 

Theorem 2.8 The ball arrangement BoptiR^ , K^) provides the densest mul- 
tiply transitive ball packing belonging to the x R space group 4q. I. 2 

{q = 2). 




Figure 6: The densest ball packing Bopt{Ri, K4) is illustrated with parts of some 
neighbouring half spheres B^, B^^ss, Bj^2t and the "base plane". The spheres 
i?A'Tg3, Bj^2t sourround and touch the sphere B^. The centres of the above 
spheres are described in Fig. 5. 
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Figure 7: The densest ball packing is described by its balls B^, B^rg^ and a part 
of the shere Bx2t . 

3 The conjecture for the densest ball arrangement 
in Thurston geometries 

The notion of the density of an arbitrary congruent geodesic ball packing in spaces 
of constant curvature is known (see Section 1). 

In Section 2 we have introduced the density function of the geodesic ball pack- 
ings generated by a discrete isometry group in a given Thurston geometry. This 
density is related to the Dirichlet-Voronoi cells generated by centres of the balls. 
For these ball packings we can formulate the following 

Conjecture 3.1 Let B be an arbitrary congruent geodesic ball packing in a Thurston 
geometry X, where B generated by a discrete isometry group of X. The above 
determined ball arrangement Bopt{Ri, K4) with density 5(i?4, K4) ^ 0.87757183 
provides the densest congruent geodesic ball packing in the Thurston geometries. 

Remark 3.2 So far the shape of Dirichlet-Voronoi cells moreover the equidistant 
surfaces of two points are unknown in some spaces (see ^T5\l , 4731/ . 4771/). 

The general definition of the density of congruent geodesic ball packings in the 
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Thurston geometries is not settled yet, but by our investigation for any "good" 
density definition can be formulated the next 

Conjecture 3.3 The densest congruent geodesic ball packing in the Thurston ge- 
ometries is realized by the above ball arrangement Bopt{R4, K4) with density 
5{RA,Ki) ^ 0.87757183. 

We are working in similar problems in SL2R and Sol space, too. 

In this paper we have mentioned only some problems in discrete geometry of 
Thurston spaces, but we hope that from these it can be seen that our projective 
method suits to study and solve similar problems ( [|22]| . [|2T]| . [|26l ). 
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